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DEFINITION OF SYMBOLS 

Latin Symbols 

Definition 

A 

algebraic substitution used to simplify an expression 

C 

integration constant used in the DuBois-Reymond equation 
when subscripted 

C 

fixed discontinuity in a state variable when unsubsc ripted 

D 

specified relationship between variables at a point of 
discontinuity 

E 

Weierstrass excess function 

F 

fundamental function appearing in the formulation of the 
multiplier rule 

g 

a function of the end points to be extremized 

J 

the numerical value of g 

l 

Lagrange multipliers associated with the end point 
constraints 

R U 

a region of 2n+l dimensional space in which the constraints 
are assumed to have three orders of continuous deriva- 
tives 

Rl 

the union of all R regions 

Ra 

a region such that the m(n+l) dimensional point 
[ tj, . . . ft *x(t 1 ) , . . . » x (t m ) ] are interior to it 

t 

the independent variable 

X. 

1 

any of the set of original state variables 

X 

the entire set of x. 's 
1 

X 

used only as X: a set not identical to x in the Weierstrass 
test 


iv 
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DEFINITION OF SYMBOLS (Cont'd) 

Latin Symbols Definition 

z. variables employed after the transformation which include 

1 both the initial state variables and the value of the inde- 

pendent variable at which discontinuities occur 

z the entire set of z.'s 

1 

Z used only as Z: a set not identical to z in the Weierstrass 

test 


Greek Symbols 
A 


e 



X 


IT. 

1 

T 

<P 

$ 


indicates the magnitude of discontinuity in the i state 
variable at the i^ 1 time point as in Ax. . 

a small number 

a Lagrange multiplier 

the entire set of \.'s 

l 

a set of numbers used in the Clebsch condition 
independent variable in the transformed problem 
a dynamical constraint 
an end point constraint 


Subscripts 

f indicates a final point 

h dummy summation variable: range 1, .... n 

i index used on a member of the set {x lt . . . , x n ) 

j index used on a member of the set {tj, . . . ,t } 

k index used on a member of the set {z^ z^} 

or dummy summation index 

v 
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DEFINITION OF SYMBOLS (Cont f d) 


Subscripts 

m 

n 

N 

P 

P 

q 

Q 

r 

s 

a 

IS 

v 

y 

0 

1 

2 

Superscript 


Definition 

the number of values of the independent variable at which 
discontinuities occur in the state variables 

the number of state variable in the original problem 

the number of state variables in the transformed 
problem 

the number of end point constraints 

the number of end point constraints in the transformed 
problem 

the number of constraints in the original problem 

the number of constraints in the transformed problem 

dummy summation index: r = i, . . . , m 

dummy summation index: s = 1, . . . , N 

index on constraints in transformed problem: 
ot — 1, . . ♦ ,Q 

index on constraints in original problem: (3 = i, . . . ,q 

index on end point constraints: v = 1,. . . , p (original 
problem) or v - 1, . . • , P (transformed problem) 

index on constraints in original problem: y * 1, . • . ,q 

indicates an initial point in the original problem 

indicates an initial point in the transformed problem 

indicates a final point in the transformed problem 


j 


indicates the particular arc under consideration 


DEFINITION OF SYMBOLS (Concluded) 


Other Notation 

* 


Definition 

indicates a quantity expressed in terms of the transformed 
variables 

indicates a derivative with respect to the independent 
variable of the initial problem 

indicates a derivative with respect to the independent 
variable of the transformed problem 


vii 



A REVIEW OF SOME EXISTING LITERATURE CONCERNING 
DISCONTINUOUS STATE VAR IABLES IN THE CALCULUS 

OF VARIATIONS 

SUMMARY 


The important problem of determining optimal trajectories for prob- 
lems with discontinuous state variables has been treated by a few authors in 
rather esoteric papers. This paper summarizes the contents of two of these 
in a fairly detailed discussion. An attempt has been made to develop ex- 
plicitly those relationships which are not contained in material dealing with 
the case of continuous state variables. 

The essential feature of the discussion is a transformation from a set 
of state variables with points of discontinuity to a set of new variables which, 
though greater in number, are continuous. The necessary conditions of 
Euler, Weierstrass, and Clebsch, along with the transversality conditions, 
are discussed in the transformed problem; the results are rewritten in the 
original variables. 

The bulk of material covered assumes that the magnitudes of the dis- 
continuities are known a priori, but the case of unknown discontinuities is 
treated in the latter portion of the paper. 


INTRODUCTION 


Most of the modern texts dealing with the subject of classical mechanics 
are oriented to those areas which form a basis for quantum mechanics. The 
recent interest in the application of classical mechanics to astronautics 
cannot be adequately treated by the texts which deal primarily with ground- 
work for quantum mechanics, and for this reason much research has been 
done into the older literature of mechanics. 

One such area which arises in the aerospace field is the necessity of 
considering discontinuous integrands in the problems formulated via calculus 
of variations techniques. This problem has been treated only very recently, 
and the formulation rests on the doctoral dissertation of C. H. Denbow which 
was published in 1937 (Ref. 1) . 



Denbow did not assume discontinuous variables but rather a problem 
situation in which unknown corner points existed. The basic idea of this 
paper was extended by Hunt (Ref. 2) to a case in which the variables were 
assumed to be discontinuous. The points at which the discontinuities appear 
are unknown, but the magnitudes of the discontinuities are assumed known. 
Mason, et. al. (Ref. 3) , have treated the problem in which both the loca- 
tions and magnitudes of the discontinuities are unknown. The work of Boyce 
(Ref. 4) extends Hunt's treatment to explicit recognition of Ihe control varia- 
bles and to the inclusion of inequality constraints in the formulation. 

In order to obtain some insight into the necessity of considering such 
problems, we note that it is often desirable to determine a so-called opti- 
mum rocket trajectory, that is the trajectory which maximizes - or mini- 
mizes - some quantity depending upon the end points of that trajectory. 

Within this framework we can treat the problem of finding the maximum pay- 
load that a vehicle can deliver to a specified orbit. 

The problem developed by Hunt covers the case in which there is a 
known vehicle, i. e. , the dry weights of the individual stages are known. By 
this assumption we know, at launch, what the discontinuity in the mass must 
be at each staging point. (Note that the thrust actually "tails off" then "builds 
up" so that thrust may be regarded as a continuous variable. Then the ac- 
celerations are discontinuous only at points where the mass is discontinuous. ) 

The problem of Mason et. al. is applicable to the case in which an 
existing rocket is not specified, but rather one in which individual stages are 
designed in such a way that a maximum payload with given liftoff thrust is 
obtained. It is well known that optimum weight ratios exist between suc- 
cessive stages; since these are coupled with the trajectory to be flown, we 
have a full calculus of variations problem. 

This paper does not add anything new to the above theories. The pur- 
pose is, rather, to expand the above work to such a level that most workers 
familiar with variational methods (at the level of Bolza's problem) can obtain 
a knowledge of the discontinuous problem with a minimum of time expendi- 
ture. All of the above papers are rather formidable and leave much of the 
developmental work to the reader. 

The notation used in Hunt's work will be retained, with the exception 
that certain conventions (such as summation on repeated indices and a sub- 
script notation for partial derivatives) wiil be replaced by more explicit 
notation for clarity. The subscript notation is probably the most difficult 
part of the study, and any simplification proves most beneficial. 
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Finally, it will prove worthwhile to carry references to a specific 
problem throughout the development in order to both illustrate and moti- 
vate the theory. The example that will be chosen is, of course, the multi- 
stage rocket problem. 


FORMULATION OF THE PROBLEM 


The problem of Mayer considers the techniques of finding a maximum 
or minimum of some function of the end points 


J = g[t 0 ,t f ,x(to) , x(t f )] 


(where t 0 ,t£ are the initial and final times; x(t 0 ) , x(tf) are some functions of 
these time points) . It is assumed that the x's must satisfy certain differen- 
tial equations and end point conditions. In our present problem we shall be 
concerned with a finite number of time points at each of which certain condi- 
tions shall be specified. ( These can be thought of as the points at which 
staging occurs, for example.) 

To illustrate the problem, we assume that we have n variables, 

Xj(t) (i = 1, . . . ,n) and m time points tj, j = 1, . . . ,m. Assuming that we 
already know the answer, as is usual in calculus of variations, the following 
diagram illustrates the x^(t) for 3 variables and 4 junctions, i. e. , n = 3, 
j = 4. 



In this case the x^(t) have been chosen to maximize a function 


J = g[ti,t2,t 3 ,t 4 ,X 1 (t 1 ),X 1 (t 2 ),X 1 (t4),...,X3(t 1 ),...,X3(t4)] 
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In general, we write 


J _ Si • • • » tyy, ) X (tj ))•••) X( tj^) ] 

where x(tj) = [ x^tj) ,x 2 (t 1 ) , . . . , x Q (t 1 ) ] etc. 

The x’s are, again as usual, subject to certain differential equation 
constraints along the trajectory. If we have q such constraints we can write 
these as 

<pp( t, x, x) = 0 (/3 = 1,... ,q) 

But if, for example, we wish to change the nature of a differential 
equation at some time point tj then we must append an additional label to ob- 
tain something like 

<p^ (t, x,x) = 0 (valid for t ^ t.) 

(Specifically, we might wish to let be a flight path equation which in- 

cludes aerodynamic effects and <p^ + * be the corresponding flight path equa- 
tion with the atmospheric terms deleted. ) It should be noted that the number 
of constraints must be less than the number of variables or there will be no 
freedom to choose an optimum path. Thus, we require 


q < n 


with the strict inequality assumed to hold. 

To obtain the range on j we relate this index to our number of end 
points. For arc number one, with tj ^ t < t^, we have j = 1; that is, all our 
equations on the first arc are labeled with a 1, and the differential equations 
are 


<i>pH t,x,x) = 0. 
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The second arc is labeled with a 2, etc. We must include the last end point, 
tjjj, for some arc, and since we have excluded the right-hand end point on 
each arc we write 


For j = 1,2 m-2, we have tj 


t < t. 


j+1 


For j = m - 1 


we have t. ^ t. . 

J J +1 


(Note that m points give m-1 arcs, so in the 
never equal to m. ) 



= 0 equations we have j 


Other constraints on our problem are the conditions which must be 
met at the points tj. (This condition allows us to specify whatever seems 
appropriate at each staging point such as initial latitude, final orbital in- 
clination, altitude of second stage cutoff, etc.) If we label our constraints 
at the tj points by we have 




• x(, m n 


= 0 


where, again, the x f s are assumed to represent the entire set (x lf . . . 9 x n ) . 
Now the range on v must be determined. This procedure essentially counts 
how many things we can fix. We may, if necessary or desirable, fix all time 
points tj, which gives m conditions. We may also fix all the x f s at t A and 
t m , giving 2n more points. Finally, we can specify all the x f s at either 
(not both) side of an intermediate time point, tj. Since we have n x's fixed 
at either side of m-2 points, n(m-2) more conditions are added. (We cannot 
fix the x T s at both sides of tj since we shall later specify the magnitude of 
the discontinuity; this procedure could yield internal inconsistencies. ) We 
obtain a total of 


m + 2n + n( m-2) = m( n+l) . 


This is the maximum number of constraints we may have. Thus, 


v = i,. . . ,p; P ^ m(n+l) . 
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These remarks preface the following statement of our problem. 


Find in a class of admissible* arcs 


Xj(t) (i — !»•••> n; tj — t £ t^) 


( 1 ) 


satisfying differential equations of the form 


( t, x, x) = 0 


(/3 = 1 q < n) (2) 


(for j = 1, . . . , m-2 we have t. st< t. + ^ 


for j = m - 1 


we have t. s t ^ t. . ) 

J J+i 


and end conditions of the form 


ip [ t m ,x(ti) x(t m )] = 0 (v = 1 p=s m(rn-i)) (3) 


one which minimizes 


J — g[ tj, . • • , t^, x(tj ) , , • . , x(t m )] . 


(4) 


We have now isolated all discontinuities at the points h. Since we 

/ . v J 

shall apply the multiplier rule, we require that the <p J' functions are of class 

C 3 (i. e. , (pjf'eC 3 ) . This requirement is assumed to hold in some region of 
a 2n+l dimensional space (n x's, n x's and one t) which we label as R^j (the 
j is discussed below). Furthermore, the multiplier rule requires that the 

equations are independent. Since this condition must hold on each arc, 
we may state independence by requiring that the matrix 


* The term ’’admissible” requires some explanation and will be defined on 
page 7 , paragraph 3. 
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ff£| 

9x. 

1 1 

be of rank q. (Remember we have q equations.) Thus, we are assured 
that we can uniquely solve for the q x.'s. 

Another requirement which must be assumed for the following theory 
to be valid is that our 2n+l dimensional points are all interior to the region 


We require an independence relationship among the ip u 's for similar 
reasons. That is, we require that the matrix 


% \ 

at. ’ 9x.(t.) 

3 13 


(5) 


be of rank p (we had p end conditions) . Also, let R 2 be a region such that 

the m(rH-i) dimensional points [tj,... , t^.xft*) x(t m ) ] are interior to 

R 2 . In this region we require that all ip v ’s and g are of class C 3 . 

To return to the subscript which was placed on R t , the raison d'etre 
for j is that we considered each individual arc on which the differential 
equations held as separate segments. If we now take the union of all R^j's, 
we obtain the total region Rj. In set notation 


Ri = 


m-1 

U R 

3=i 


ij 


(m tj's give m-i arcs). An arc in Rj can be defined as an arc x^(t) over 
the range t A ^ t ^ t m which is continuous except possibly at a finite number 
of points, t. (j = 1, . . . , m) ; x^(t) consists of a finite number of arcs, and 
x^(t) is assumed to have continuous first order derivatives (x^(t)eC 1 ). 

The term "admissible" arc can now be defined. 

Definition: An admissible arc for our problem is a set [ t* t ] 

and an arc interior to R lf all of whose ends and corners 
[ t. , x( t.) ] are interior to R 2 . 
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By this definition an admissible arc gives well defined values to 0 , g, 

ip v , and J. 

This information is preliminary. We turn now to the transformation 
of our stated problem into a form tractable to classical methods. 


TRANSFORMATION OF THE PROBLEM 


As is well known, the theory of variational techniques cannot handle 
problems involving discontinuous variables. These could be applied if a 
technique existed to make all variables continuous. The observation of 
the diagram on page 3 and a little reflection might give rise to the following 
concept. The variable Xj is obviously discontinous at the point t 2 . But if 
we were to regard Xj on the first arc and x A on the second arc as two 
variables instead of one discontinuous variable, then we have a problem in- 
volving only continuous variables. We bought this condition by adding 
another variable, but that is not a strong penalty since we already know how 
to treat the multiple variable problem. 

The discontinuity of Xj at still exists, of course, since we have not 
altered the physics of the situation. We must eventually find a method of 
connecting our path across this point, but, as will be later demonstrated, 
the natural corner conditions will do just this. 

This, in essence, is the philosophy behind the manipulations that 
follow: transfer from a problem involving discontinuous variables to one in- 
volving an expanded set of continuous variables. 

We initiate the procedure by transforming the time. To do this, in- 
troduce a ficticious time r and define t(-r) as 


t = t + (ti +1 - ti) r (j = 1 m-1) 


with t having a range 0 ^ r ^ 1. For our first arc we have 


ti ^ t< t2, 

so t =tj + (tj - ti)r which gives j = 1. For our last arc, 


t 

m- 


1 


£ t < t 


m’ 
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SO t = t 4 + (t - t )T, 

m-1 m m-1 


so j = m - 1. 


The variables x. are functions of t, and we introduce new variables 
and call them z's, whioh are functions of r. They are defined in terms of 
the x's by 


z 


(j-l)n+i 


(t) = x^t) . 


( 6 ) 


That subscript on z looks a bit odd, but the original approach shows that it 
can be explained. On the first arc (tj ^ t == we have variables x*, . . . , x Q . 
For this arc, as noted above, j = 1. So on the first arc 


z^r) = V 1 ^ (i = 1, . . . ,n) 


Now we must introduce new z’s for the second arc, and since we have used 
up Zj z n we start the second arc variables as 


Z n+i (r) =x i (t) * 

Therefore, 

z .(t) = x.(t) (secondare), 
n+i i 

which can also be written as 


Z (2-l)„ + i (T ' 


Since j = 2 on the second arc 


Z <j-l)n + i (T) =x i (t) 


, we have 


(second arc) . 
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Continuing in this way we find that 


Z (j-l)n + i <T)=X i (t) arc >- 


As noted above, for the last arc, j = m - 1. So we obtain our j range as 


j = 1, .. . , m-i. 


To see how many variables we really have introduced, we set j and i to 
their maximum values and find 


Z (m-2)n+n Z (m-l)n’ 


which tallies with the figure where we had n = 3, m = 4 giving z 9 . A simple 
count on this diagram shows nine variables. (Here the first cloud appears 
on the horizon; for it is not unusual to use 12 x's and 5 staging points in a 
trajectory - which means 48 variables to consider. ) 

Now we consider our discontinuity in terms of the z's. Our variables 
are defined at the left end of an arc, which is the right side of tj denoted by 
The value - which is unspecified - of x^(t) at the left side of tj (i. e. , 
right end of the previous arc) we shall call t. . Here we consider that 

\ 

t. + = t. + e 
J J [ 

> e small 

t.“ = t. - e 
J J 

J 

Directly from the definition of the z’s (eq. 6) we have 



x.(t. + ) = z 


(j-l)n+i 


( 0 ). 


At the other end of this interval we have the same x^ evaluated at tj+i . We 

s till have z, . JV . since we don't start our new set of z's until we move to 
(j-i)n+i 

the right of t So 
j+i 
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x.(t. .“) = z. . . . .(i) . 

1 j+1 (j-i)rH-l' 

For the value of Xj at the right of t i. e. , Xj(t + ^ + ) , we note that we are 
into the next set of z's; thus we add n to our previous value giving 


x.(t. + ) = z . . . .(0) =z. .(0), 

1 j+1 (j-l)n+n+i jn+i 


Adding the last two equations we have 


z.. .\l)+x.(t. , + )=x.(t. ") + z. .(0) 

(j-l)n+i i v j+1 i' j+1 jn+i v 


or 


z . . .(1) + [x.(t. . + ) - x.(t. , ) ] = z. .(0) 

( j-l)n+i v l j+1 l j+1 jn+i 


= z U-l)«»i (0K 


(7) 


The term in parentheses is our known discontinuity of the at the point 

Vi* 

Now this is good for any x, so i = 1, . . . ,n. Our last discontinuity in 
the x's occurs at t 
must have 


m-1’ 


which in our equations shows up as 




So we 


j+l=m-l— - j = m - 2. 


Our first discontinuity occurs at 1^, so x. (t^) ► x. (t. .) — j = 1. Thus our 

range of j in equation (7) is J 


j = 1, .. . ,m-2. 
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This defines our first (m-1) n variables. We now add to this the time 
points ti and make them variables of our problem. Since we have numbered 

• l . x l * xl_ „ _ .1^1 


our z's from z lt 


'(m-l)n 


we start numbering the new z variables (which 


are really just the tj's) with z ^ • • • 

Then ^ 


, z. (we have m t.'s) . 

’ ( m- 1) n+m j 


t. 

J 


Z (m-i)n+j* 


Since our tj's are constant we have 


' ( m- 1) n+j 


= 0 (j = 1 m) 


( 8 ) 


where z' = — for any z. These are all the variables that we shall need, 
dr 

We have collected a total of 


(m-1) n+m = N 
variables. 

The function g, which we wish to extremize is 


J — Si tjj •••! » • • • » X( t^) ] • 


The values of all the t.'s in this equation and the values of Xj's at the corner 
points (at one side or ihe other of the corner points - and we chose the right 
side for every point except the last one) have now been expressed in terms 
of the z^(0) , z^l) , which we shall denote by z(0) , z(l). Thus J becomes 

J* = gs=[ z(0) ,z(l)] 


where the * indicates a functional dependence of z. 


Next we transform the (f>^ *s. As in the case of the z's there will 
( i) ' 

be many more 4>\ 's in the representation in terms of t,z, z 1 than there 


were in terms 


i p o a in 

Of t,X, X. 


The new functions are defined as 




(9) 


The justification for the ( j — 1) q+/3 subscript is much the same as before. On 
the first arc, j = i; we range /3 from 1 to q (q constraints). On the next 


arc we start by numbering <p 


* = (P* 


q+i ^ (2-1) q+i r ( j-i) q+1 


= ft 


(j = 2 on the second 


arc) , etc. Again, for m time points tj we obtain m-i arcs. So in equation 
(9) we have 


j = 1, . . . , m-1 

js = i,...,q 


and we still have 


t. + 
J 


VV T - 


This condition accounts for our original differential equation constraints. 
But equation (8) also produces additional constraints of certain z' 's, so we 
can write 


<P 


* 

(m-i) q-t-j 


(t , z, z') 


z 


(m-i) n+j* 


(10) 


The subscript of the left side was cnosen to pick up directly after the last sub- 
script of equation (9). That is, we insert our maximum value of j(= m-i) and 
/?(= q) in this equation giving 


» 3$ C j 

9 (m-i-i)q+q 9 (m-i)q* 
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and start equation (10) at jjq+j* 


The right side subscript comes directly 


from equation ( 8) . Note that the two subscripts do not agree numerically. The 
only variable subscript in this equation, j, to account for all the time points tj, 
is 


j !)»••) m. 


We had previously defined (m-l)n+m = N variables and now we have 
defined a total of (m-1) q+m constraints. If we subscript all of our <£* 's with 
ana, we see that the first (m-1) n of these 


0* (a - 1,. . . , (m-1) q) 


are defined for 0 ^ r < 1 (excluding the right end of our interval again) . The 
second set 


<p*x (a - (m-1) q+1, . . . , (m— 1) q+m) 


are defined for 0 < t < 1. (It does not matter which value of t is plugged 
in to equation ( 8) . ) The total number of constraints we label as 


(m-1) q+m = Q. 


For the end conditions we had ip^ = 0, where v = 1, . . . ,p^ m(n+l) . 

We can now get - i.e. , transform ip v into dependence upon z, t variables - 

just as we got g* from g. This gives us our first p conditions. But there 
are others. 

In equation (7) we define 


x i<W ' x i‘W 


= Ax.., 
ij 


(Note that our discontinuity here is at the right end of the j n arc. Thus it 
occurs at the j+l’st point t, + ^. This accounts for the difference of the j 

subscript across the equation.) Equation (7) may now be written 


if. .. . = z,. L .(i) +Ax„ -z .(0) = 0. (11) 

r (j-l)n+i+p (j-i)in-i' i] (j-l)n+n+i v 


The left subscript starts with j = 1, i = 1, so we pick up at as we 

should since above we used up pip's. For the last point we have tj = t m _^ 
so our last arc is the j-2'nd one giving a range to j of 


j = 1, . . . , m-2. 


i, as usual, refers to our variables x lf . . . , ^ so we have 


l 1 1 • • • , n • 


Putting j and i to their maximum values we find we have accumulated 
(m-2-1) n+n+p = (m-2)n+p. 

Our final two end point constraints are on the end values of t. 

Since our last constraint was 26* we have 

r ( m-2) n+p 


^(m-2)n+p+l Tl ^ 

^(m-2) n+pH-2 T<i * 


Thus our new (transformed) end conditions are 


ip* ,t 2 ,z(t 2 )] * 0 


where 


v = 1, . . . , (m-2) n+pt-2 = P. 


The original problem has now been transformed into the problem of 
finding, in a class of arcs defined by Zj(t) satisfying <£* ^ o and ip* = 0, 
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Illllllll 


an arc which minimizes J* = g^[ z(r 1 ) , z(t 2 ) ] . This now satisfies the 
hypothesis of Bolza's problem, and we can apply the standard formulation 
as given by Bliss (Ref. 5) . The proof that the solution to our original prob- 
lem has a solution which is equivalent to that of the transformed problem when 
written in the original coordinates is not obvious from a mathematical view- 
point. The argument which justifies this statement is given in Reference 2, 
and it will not be reproduced at this time since no real gain would come of 
it. We shall now proceed to examine in detail the derivation which pro- 
duces the Euler equations, transversality conditions, Weierstrass condi- 
tion, and Clebsch condition for the discontinuous state variable problem. 


THE EULER EQUATIONS AND TRANSVERSALITY 
CONDITIONS 


Since the problem of Bolza is now applicable to our task, we can 
apply the multiplier rule as given by Bliss (Ref. 4) . An admissible solu- 
tion, E* , of the equations 0* = 0 ( a= 1, . . . ,Q) defined on [ Tj, t 2 ] , is said 
to satisfy the multiplier rule if there exist constants A.q , (v = 1, . . . ,P) 

not all zero and a function F : ‘ ( t, z , z ' , A* ) = ^ A* <p* with multipliers 

Q!=l 

A*,(r) continuous on [ tj, t 2 ] except possibly at values of r defining corners 
of E* where they have well defined right and left limits, such that the 
equations 


(k= 1, . . . , N) 

are satisfied along E* , and such that the equation 

2 

, +*§ dg'= 


P 





N 

- 2 

k=l 


, 9F* 

9z ' 
k 


N 


k=l 


dI + 2 fr dz 


9F* 

dz! 


f 

J 9z, 
Ti k 


dr + C, 


( 12 ) 


(13) 


holds at the ends of E* for every choice of differentials dr 1 ,dzj(r 1 ) ,dr 2 , 
dzj ( r 2 ) . Every minimizing arc, E* , for the pro blem of Bolza m ust s atisfy 
the multiplier rule. 
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The function F= defined in the above statement can be written by use 
of equations (9) and (10) as 


f: 

j 


= t 

/3=1 


A* 

( 3 - 


i)q + /3 (T) ^(3-l)qV T ’ Z * Z,) 


(14) 


+ 


m 

y ^ 2, * 

. (m-l)qfr ( m- 1) n+r 
r=l 


where we have split our constraints into the two groups defined by equations 

(9) and (10). The j subscript affixed to P* is (mechanically) necessary 

because of the free j subscript on the right hand side of the equation; its 

physical significance corresponds to the possibility of F* differing from 

pi; by variations in the <b* f s. Note that j is not summed, (Obtaining a 
j+1 

non-subscripted j would require an additional summation.) The equations 
produced from application of equation (12) to the function F* have different 
characteristics depending upon whether we consider the first (m-1) n of 
these equations or the last m equations. 


Application of the Euler Equations to the 
First (m-1) n Variables 

For the first (m-1) n equations, we begin by consideration of the 


term 


8F? 81* n 

1 _ 3 V 9F 


.1 _ V JL? 

«•» t-l P)v 


dx 


8z' 8z' . . u ,dk dz* . 

k (j-l)nfi s=l s (j-i)n+i 


But 


x i (t) = z (j-l)wi (T) 


so that 


k i (t) = 2 (j-l).widt- 
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Since 


t=t i + VV T - 


we have 


Then 


giving 


dt 

— = t - t . 
dr j+i j 


. . = (t. - t.) x. 

(j-l)n+i j+l j 1 


dx. 


dz' ... . t. . - t. 

(j-1) n+i j+i j 


In the above sum we obtain a non-zero term for only one value of s. 


9F: 

L 


9F 1 

9x. t. - t. * 
i J+i J 


From equation (14), by direct differentiation, 


9F ; ; q 
1 = V 

8k i / 3=1 8k i 


"hum 


and since 


Fj[T(t) ,z(x) ,z'(x,x)] = F[ t,x, x], 


Then 


18 



we have 


3F* Q 

3z' = 2, 

k j8=i 


(j-i>W 

t. . - t. 
3+1 3 




(i-Dq+g 

ax. 


The next term of equation (12) which we develop is 


8F* 


n 


- 2 

k=l 


3F* dx 
j_ s 

8x dz. 
s k 


8F* 



(k = 1, . . . , (m-l)n) 


dx. 

(again, x depends only upon z and 

S Q Z . . . v 

(j-l)n+i 


1) 


aE* 


ax. 

i 


q 

= 2 

0=1 


dp 


x* 


(i-i)q+/3 


(j-i) q+/3 


ax. 

i 


Equation ( 12) thus reads 


q X* 

j (j-i) q+£ 

J3=l Vl " *3 


Tu-Pttfl. 

ax. 


" pk X * (i - 1)q+ P 


9 ^li-l)q+/3 

ax. 


dr 


(15) 


C (j-i) rri-i 


= 0. 


We now transform multipliers to convert equation (15) into a more 
useful form. Since X$ is a constant, we have the definition 


X 0 - Xi<, 


(16a) 
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and the transforms of the multipliers A* (r) are defined by 



(t) 


( t - t. 


X *1-W T > 

Vi-*j 


(j-i)q+/3 


It. . - t. , 

J± i . -Jj 


t. . - t. 
3+1 3 


(16b) 


Note that in the right hand term we have simply replaced r by 


that is, A* ,, LO 
(3-l)q+j3 


t - t. 
L 


t - t. 

1_ 

t. . - t. 
3+1 3 


it. . - t. 

V 3+1 3 1 


indicates a functional dependence . Note that 


t , - t. 
3+1 3 

given j. 


is a constant multiplying all multipliers by the same value for a 


Equation (15) now reads 


k e <*i » Sx i "- lhl+i 

3 


where we have substituted 


(17) 


<P 


(j) 

P 


= <P 


(3-D q +P ' 


To explain the change of the integration limits, we note first that since 


t. 

3 




we have 


t| =t| n =t. 

Ti T=0 J 
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which explains the lower limit. The upper limit, r , is 


r 


t- 




which is actually our upper limit. But since the integral is taken to any 
value of the independent variable, now t, we simply set 


t - t. 

1 1. 


t. , - t. 
3+1 3 


We can now write our function F as F(t, x,x,\) with 


F - t v> 


/ 3»1 


P [3 


where it is now understood that we restrict ourselves to the arc 

t. < t ^ t. . so that 
3 3+1 


= ^ ’ A /3 = X i3 j) f ° r Ij ~ 1 < tj+i 0 = I,- • • .m-1). 


The range on j, as usual, is due to the m-1 arcs generated by m points 
t*, . . . , t^. Equation (17) now assumes the fairly simple form 


§r* - f ~~ dt - C . = 0. (t.^t<t. ). 

ax. j 9x. (j-l)n+i J J+1 

3 


(18) 
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Application of the Euler Equations to the 
Last m Variables 


Thus far we have applied equation (12) (a portion of the multiplier 

rule) to the set of variables corresponding to our original x's. But we have 

a number of other variables to consider; these are the points t ls . . . , t^, 

and we now investigate the application of the Euler equations to these. This 

will result in determination of the multipliers A* . . 

(m-1) q+r 

To reiterate the applicable equations, we have, first, 


K 



8F* 


9z, 



(k = 1, . . . , N) 


where N (= (m-l)n+m) is the total number of constraints. In deriving equa- 
tion ( 18) we used up only the first (m-1) n variables (i. e. , the original set 
of x's). To proceed to the last m equations, we use the DuBois- Reymond 
equation ( 12) (which is really a generalized form of the Euler equations) . 

As before, 


Ft 

J 


q 

2 A* d>* + 

4 (j-i)q+/^(j-Dq+/3 


m 

(m-i)q+r Z (m-l)n+r * 


The first derivative yields 


3 Ft 

9z (m-l)n+l 


= A 


* 

(m-1) q+i ‘ 


Note that a specific value of r has been chosen - namely 1. We shall ex- 
amine the first of these new equations in detail. 


3 Fj 

The next derivative to be taken, , for k = (m-l)n+l, requires 

k 

more work. Recalling that z ^ m _^ n+ ^ is just another label for t* we have 
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9Fr 


9F* 


9z (m-l)rH-l 8t > 


' wT \ W T) *0-« w> ( tW ' z(x ' tK: 


By definition of t 


'-V VV T ’ 


and the only place we have tj is when j = i. In that case 


t = tj + (tg - tj) t . 


Then if we set i = 1, P* becomes 

J 


9F* 


9t< 9t 


-& {IKK = l x 


£1 » f* 


% 

8t » 


l*, 

>3=1 P 


d<p* n 9tf>* 9x. n 9tf>* 9x. 

_£ ai . y _£ _±_ + y ^ 

9t 9ti /-< 9x. 9t t 9x. 9tj 

^ i=l 1 1 i=l 1 1 


From the equation given for t with j = 1 we have 


9t_ 

9t< 


= 1 - r . 


'(k)) 



I 


I II III I 


II II 


III 


Also, since the x.*s and t's are independent variables we have 


9x. 

i_ 

at. 


= o. 


Then 


^ = T \* 

at, U A « 


at 


/3-1 


% n a** ak. 


at 


As commented previously 


x. 

1 


dr 


z 


i dt 


and 

T = ~~ t (for j = i) , 

h ~ h 

so 

dT = 1 

dt tj - tj 

and 


8x. 

i_ 

at. 


a 

atj 


k ~ h J ( *2 — t A ) ' 


Finally, 


8F» 

at A 


arx= 


q 

- = y \* 

8z (m-l)n+l /3=i ^ 




i=i dk i 
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9 pc 9 JiVr^ 

Inserting the expressions thus derived for — and in the 

DuBois-Reymond equation, we find (m-l)n+i 1 


A* 


( m- 1) q+1 


t q 

Ti /3=1 


90* 

—1 

at 


(i - r) 


n a<p* 

+ y -£ 

Li a> 


i=l 8X i 


dr - C 


(m-1) n+1 


(19) 
= 0. 


This equation, unlike those above, does not restrict the minimizing 
arc and merely determines the multiplier i)q+i^ T ^ w ^ ic ^ I s associated 

with the equation z' . . . = 0. The usefulness of this derivation is not 

(m-1) n+1 

immediately apparent but equation ( 19) will be of value in later work. 


First Use of the Transversality Condition: 
Application to the Coefficients of dz( m -i) n +j 

The transversality condition, equation (13), must be written in ex- 
panded notation. Thus, 


dg* [z(t 1 ) , z(t 2 ) ] 



9g* 

az s ( Tl ) 


dz s (Ti) + 


N 

y 

< 9z 

s=l s(t 2 ) 


dz s (T 2 ) 


where the sum is taken over all z variables, s = 1, . . . ,N(=(m-l) n+m) . 
Furthermore 


N df 

d^[Ti,z( Tl ),T 2 ,z(T 2 )] = Yj 9z " ( Ti ) dz s f T l> 

s = 1 s 


N 8#= 

+ Z p 


dip' dlfr 

^8z„(t 2 ) dZ B (l * ,+ ^T dT « + -^ dT *- 


s=l s' 
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Since transversality must hold irrespectively of the values chosen 
for the differentials, the individual coefficients must vanish. From equation 
(13) we gather together the individual pieces of the dz^ coefficients and 
write expressions at the end points as 



T 2 


dZ k (T 2 ) 



dip* 

r _V 

az k (T 2 ) 



T 1 


+ x 0 


9g* 

9 z k (Ti) 


p dip* 

+ y i _ l _ 


= o . 


Summations on s do not appear since we chose the specific coefficient s = k. 

Furthermore, the sum on v extends over all end conditions 1 P(=(m 

-2) n+p+2) . The above equations are decoupled at the points t*, t 2 since 
dz k ( T i) and dz k (T 2 ) are independent insofar as transversality considerations 
are concerned. 


We now restrict our discussion by deciding on a value of k to be in- 
vestigated. As was done in the derivation of equation ( 19) , we use k = 
(m-i)q+l. The first specific equation in this derivation was 


8 Fi' 

L 


8z 


= \* 


( m-i) n+i 


(m-i) q+1 * 


The above two transversality conditions now take the form 


V (m-l), + l (T 2 ) + 4 te (m _ 1)q+1 (T 2 > 


Z 1 , 


dip* 

Y v 


dz 


= 0 


v=l (m-i) q+1 (t 2 ) 


(20a) 



X (xn-i)q+l (Tl) +X ° Bz 4 ) 


(m-l)q+l 


P 

+ Jj 


dip* 

V 


v=i 9Z (m-i)q+i^ Tl ^ 


= 0 . 


(20a) 


( Note that our variable z . cannot enter ip* . and 

( m-1) q+i T ( m-2) n+p+1 

^(m-2) n+p+2 so that our sum could run onl y to P -2 » but it makes no dif- 
ference if we carry it to P.) 

We now choose A.^ and £ p in such a way that 


*5 


3*£L 


9Z (m-l)q+l^ r ^ v=i 


P 

+ 2 l v 




9z (m-l) q+1 ^ 


= 0 


so that our second transversality at t* gives 


'1m-l)q + l (T ‘ ,=0 - < 20b > 

Returning to equation (19) we can choose our upper limit of integration to be 
so that the integral vanishes. This leaves 


^(m-l)q+l^ Tl ^ C (m-l)n+i °* 
but since this \ at T\ vanishes we have 


C, = 0. 

( m-1) m-1 


27 



Equation ( 19) then becomes 


A* 


(m- 1) q+i 


T q 

M-f Z# 

r, /;=i p 


% n d<j> 


jti 9k i ^ - *!>' 


dT = 0. 


If we now choose our upper integration limit as r 2 , we obtain an ex- 
pression for A* .v . (to) in terms of the integral. But our first trans- 
(m-l)n+l 

versality condition gave another expression for A* m i) n +i^ T 2^ and we can 
eliminate this Lagrange multiplier between them giving 


An 


9g» 


9z 


(m-i) n+1 


(t 2 ) 


P 

+ 2 A- 


9# 

r i> 


9z 


y=i (m-l)nfl 


(t 2 ) 


t 2 q 

*/ 2$ 

r i /3=1 


90* n 90* 

_£ ( i _ T) + V — 2 . 

at { } L 


9k i 0-2 - * 1 ^ 


dr = 0. 


( 21 ) 


Equation (21) is essentially a modification of equation (19). The in- 
formation derived by the first transversality condition has been incorporated 
into equation ( 19) to obtain (21). In order to obtain directly useful information, 
further development of equation (21) is required. This development is initiated 
by evaluating the integral in this equation in the next derivation. 

To initiate this third segment, we note that by definition of t, 


_d_ 

dt 


1 _d_ 

h. ~ h dT 


(for j = 1) , 


which gives 


x. = 

l 


1 

k - k 


z! 

i 
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and 


x. = 


i ( tg tj) “ i 


Previously we required that our minimizing arc, E* , had no corners 
between any pair of the points (t , tj ; thus we can differentiate the follow- 
ing expression (which is written simply as an inspired guess) . 


&< 2 


|8=i 


A* <t> * n 9d>* z! 

J-L _ y x * Z£ i__ 

<2-ti iti' P 


. I 1 r n 80* z' 

= iiy X * A * _ y _£ — — 

dt V=ip *' 3 i=i ? 3k i •» - *»j 


where we have converted from t to t via absorption of ^ - t 1# 


The right hand side is now expanded - by the chain rule, and we pick 
up a second summation on the first terms since we apply the chain rule 
through the x variables: 


\ 


A V 

dt V=i 


n 8 0* 


^ ^ 


i-i n 30* 


n 8 0* 


9 <P* 0 




80* 
— $- 
P 9x 


- A: 




'j8 3x. (ta - t A ) ■ 


q n I 80* 80* 

>= E £ < X* 

j8=l i=l 


'/ 3 9x. i jS 9x. i 

l l 


/ 




(above equation concluded on following page) 
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+ A* 


90* 


+ K<P*o ~ K 


8 cb* 

— ( — @. 


p at T /3 r /3 i dt p/3 ak. I 


d<p* 

_ A - . p \ 
X P ax. i| 


9. n 

■L I { 

0= 1 i=l 




3 


x + A* — ^ + A* <f>* - x — 

A /3 ax. i A /3 at A /T /3 i dt 


d Jjl 

A* — ^ 

P 8x, 


where the relationship between z\' and x^ derived above was used, 
ferentiation of equation ( 15) , wi^h j = 1, we have 


_d_ 

dT 


j. 

2 


3 


3 


9x i 


y .. 3 

0=1 * Sx i 


or 


_d_ 

dt 


q a <p* 

li. 

_/3=l 1. 


- 2 


d<p * 

v* Zl 

/3=1 ? 9X i 


Insertion of this expression into our above equivalence gives 


_d_ 
d r 



n d<p* z[ 

y A , : ZA i 

i=l ^ 9k i (t 2 " fc l> S 




d<t>* 

ZA 

at 


X>\< (b* 


A. 

2 

0=1 


3 

V vl> — f— 

Vi' 

/3 at 


where the last step follows, since <M vanishes along an extremal, 
set 


q 

a- 2 

0=1 


. V w 3 _!L_ 

t 2 - tj £ ax. (ta - t t ) 2 


for algebraic convenience. 


By dif- 


We now 
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Then 


^ y x * 

dr L a 


0=1 ? 


at ' 


We are now in a position to begin evaluation of the integral contained 
in equation (21) . The last equation in (21) gives 




9g* 


+ Ij £ i 


8( » 


8z <m-l)»i <T » ) ' 2? 8z <m-l)n + l (T * ) 


/ 

Tl 


z\ 

1 


, . q n 9d>* 


dr = 0. 


The first portion of this integral can be written as 


7 2 


/ ^•(l-r)d T = A 


Tl 


T 2 t 2 


- / TdA • 


Tl Tj 


Setting u = r, dv = dA, we obtain via integration by parts 


J2 




Tl 


- tA 


t 2 T 2 

+ | Adr . 
Tl Ti 


The total expression for equation (21) is now 


9g* 


9z 


(m-1) n+1 


(t 2 ) 


P 

+ 2 


dip* 


dz 


v=l (m-1) n+1 


(t 2 ) 


+ [A(l - t)] 


t 2 t 2 

+ f AdT 
i Tl T t 


t 2 I q n dtp* z\ 

+ / < Z Ik — L 


Tj Uti “i <% - V 


'dr = 0 . 
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Substituting the value of A back in the integral where it appears gives 




_9 g* 


P 

+ Yj^-v 


dip* 

v 


° 8z (m-l)»l (T ^ ' £T 8z <m-l)n + l (T >> 


+ [A(i - t)] 


t 2 


To <1 \* 0 * To n q 30 * z! 

♦/ Z +i d -J 

Tj /3=1 2 1 Ti 1=1 /3=1 ^ 1 w 1 


8g* 


t 2 n q 30* z! 

+ / 2 2^ if Tt, Ji,)* dT = x? 5T .. <t 2 ) 

Ti 1=1 /3=1 p i w 1 (m-1) n+1 ^ 


P 

- 21 k 


3 ib* 

v 


,=1 8z (m-l)n + l (T ’’ 


+ [ A ( 1 - T)] 


T2 
T 1 


+ 


/ ! (iFy 


dT = 0 . 


Again using the fact that our 0* ’s are satisfied along an extremal, i.e. , 
0* = 0, this becomes 




8g* 


P 3 ip # 

Zv ” 


3z, j ( T 2 ) < P 3z , (r 2 ) 

(m-l)n+l 11 v=i (m-l)n+l z 


+ [A( 1 -t)] 


t 2 


= 0 . 


Inserting A (as was done above during the integration) gives 


q n 

z z 

/3=i i=l 


rt 


Ar’,= (h 5'.' 


d(b* z . 

v . _£ L 

A-f, , , 2 


V'Pik-W p (h-W 


(1 - t) 


T 2 


+ *8 


3g* 


ZA 


3d)* 

V 
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9z ^ ,,(t 2 ) ^ v 9z / . j (r 2 ) 

(m-1) n+1 1 v=l (m-1) n+1 z 


= 0 


L 



The only portion of this equation evaluated at Tj gives rise to 


q n 

-l l 

j3= 1 i=l 


(h * 

n 


- A* 


d<p* z . 

1& L 


(_ P P (h ~ ti) /3 8 x. (tg - t A ) 2 J 


= 0 


Previously we chose Aft ,& v in such a way that 


M 

A o 


(m-1) q+i 


P dip* 

+ 1 *, 


Sz, (n) " to, n„ J .,< T l> 


= 0 


(m-1) q+1 


We can add the last two equations to give 


„ , P dip 

. gg* . . + v ‘ 


Z 


(m-1) q+l ( Tl) ' 8Z (m-Dq+l (Tl) 


q n 

-Z Z 

/ 3=1 i=l 


\*<p* 

n ■ n 


- A 


8*| zj 

5.C ■ * ■ 


■jS j8 (ti - t t ) /3 ax. (tg - tj)' 


Ti 


= 0 


We can now make the identifications 


'(m-1) q+1 ** ’ 


At 


_ Aq , 


JL 


t2 " ti 


= A^ , and 


z i 


*2 "ti 


= X. 
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and obtain 


P dip q n 


.(lMD , % • 

X P % ' X p 3x, X i 


= 0 . 


We had defined for our first arc 


F - 

/3=i P P 


so the above equation becomes 


x ^-+ f 1 1 . 

M 8*1 " 


n 


9F 

9x. 

l 


= 0 . 


(22) 


This equation applies at our initial time point. The right hand end 
point, 1^, has a similar relationship (derived in the same way) which is 


P 9 4 > 

. u s . y i L 

0 9t ^ v at 

m v=l m 




9F 

9x. 


m 


= 0 . 


(23a) 


( The evaluation at upper and lower points as in the last two equations indi- 
cates whether the term is to be added or subtracted.) 

The rest of the equations come from the observation that the evalua- 
tion of our basic equation at t 2 yields 




9ffi 




P 

+ Ti l 


dip* 

V 


v=i ^ Z (m-i)n+i/ T2 ^ 


= 0, 
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since the first term vanishes identically at t 2 (t 2 = 1) . This term would 
vanish in any case since our x variables are defined on the interval 
t. < t < t. or Tj ^ r < t 2 and <£* = 0 along an extremal. Then 


A* d<p* 

<£* - A* — ^ 


q n 

^ 84 i 


t 2 


= 0 . 


Now the next equation for the point tg, i.e. , i)q+2 

equation 


q n 

-2 2 

0=1 i=l 


* 5 *> 


*3 


z! 

i 


*1 T i 


0 8X. (tg - tg) ‘ 


= 0 


Addition of the last three equations yields 


q n 

2 2 

0=1 i=i 


JjL 


<p* - 


X I 




(^ - tj) r 0 (ta-tj) 2 8x. (ta - ti)' 


t 2 


q n 

■2 2 

0=1 i=l 


1_ . _j3 

-e_ </>* _ — £ 


A* 80* Z 




h~ h p h ~ h 8x. - ta 


+ A 


8g'~ 


P 

+ z^ 


dip* 


° 8z (m-l)n+l (T2) 8z (m-l)nfl (T2) 


= 0 . 


Making the same identifications between the z representation and x repre- 
sentation, as was done above, and recalling the definition of F, we have 


P dip 


9£. + Vi —iL . 
a t- + L x v atz 


8tg 


v=i 


n 


8F 


F - Z k - 

H. l 9x. 
i=i 


<2 + 


= 0 , 
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The above equation could just as well have been derived in terms of any 
other point 7\ , so 


_ P dip 


. v 8t. 
v=i j 


- & 


9F 

9x. 

l 


= 0 (j = 2,. . . ,m-i) . (23b) 


Second Use of the Transversality Condition: 
Application to the Coefficients of &r\ and di^ 


The transversality condition, equation ( 13) , has not yet been played out. 

We derived our first set of conditions by equating to zero the coefficients of 

dz, jv .(To) and dz, .. (i=l,,..,n). Another set of condi- 

(m-i)n+i 4 (m-i)n+j 1 

tions could be derived by equating to zero the coefficients of dr* and $t 2 . 

These are associated with the end conditions ip*__ ^ A »r i = 0 and 


Y (m-2)n+p+2 
pliers i 


ji* : 

r ( m-2) n+p+1 

= r 2 - 1 = 0. Such a procedure does yield values of the multi- 


and l. 


( m-2) n+p+1 ( m-2) n+pH-2 * 

problem at any other point we really do not require them. 


but since these do not come into our 


Third Use of the iransversality Condition: 
Application to the Coefficients of dzj j_2) n +f 

Another important set of conditions, however, does come from set- 
ting to zero the coefficients of dz.. . . ,(r») and dz,. . . .(rd 

(]-l) n+i 2 ( j-l)n+i v 

( j = 1, . . . ,m-l; i = 1, . . . ,n) . We shall investigate the case of j = 1 (i. e. , 
the first arc) in detail. From equation ( 13) 


7 N 

N _ 


, , V A 

dr + l ,, dz 

_\ k=l kj 

k=l k 


+ *$<&< + l v df* = 0 . 
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L_ 


As before. 


and 


and 


which, 


Thus 


dg * = !i sfw d ^ (Tl)+ d V Tz) 


N dip* 


N dip * 


K = ^ dz = (Ii)+ 8 » b (t.) dz = (T!) 


dip* 

+ — ^ dri+ - — - dr 2 
9Tt 1 9 t 2 1 


m 


P ! 2 i X5 ( j _i) q+/3 ^j-l)q+/3 + ^(m-1) q+r Z (m-i) 


n+r 


for j = 1, gives 


m 


F j ^ ^0^/3 + ^-'.^'(m-1) q+r Z (m-l)n+r 

J /3=i r=i 


8E y x * d 3 y 

8z i " |i * 3z l ' |i ..i * 8i s 3z l 


y . ? w za *: 

# * 8i i 3z i “ £i * 3i i dt 
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For j = 1 we have 


T 


t - tj 

k-k 


dr _ 1 

dt _ % - t 4 ; 


then 


9F* 

8z! 


q 


l 

/s-l 


dd>* 

-A 

dk. 

1 


1 

k-k ‘ 


Transversality now reads 


q d<p* . ' 


+ 4 


8g* 

8z.(r 1 ) 


+ 4 


8g» 

8z.(t 2 ) 


+ 



9 ih* P 9 it* 

—E- + y , —v- 

8z i (T i ) £l az i (T 2 } 


Now equation (11) with j = 1 gives 


0* = z.(r 2 ) + Ax - z (r A ) = 0. 

l+p i * li n+i 1 


In our summation on tp* above, the last two terms (i. e. , Tj = 

ib* . = 0 and r 2 - 1 = it. . = 0) contribute nothing. If 

r (m-2)n+p+l i r (m-2) n+p+2 

we write our summations as 
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P 3 ip* P-2 dip* 

y i jjl = y i . v 

Li v n z 4J v 


p_2 9^J* 


3z "V 9z. ^ 9z. ' ^ . V 9z. 

i*=l i t>=l 1 v=l l K=p+i 1 


V 9z. + 


2 or 


then the last sum yields simply 


P-2 9; i* 

y * _ 2 L_ = i 

4' v 9 z.(t 2 ) i+p 
y=p+i i' z 


Note that we do not obtain the multiplier i . here since our derivative is 

n+i+p 

with respect to z^. Thus, transversality now gives 


q ^ _j. 


at, t, - 1, 


. x * - 3g* - . x * 9 gf . . 

0 9z (r 2 ) 0 az (Ti) 


p dip* P dip* 

+ kl V 8 V T 2' + ^ 1 ’’ V’*' + ii+ P = 0 


(remember the last term is associated with r 2 although it is a constant) 
giving rise to the two equations 


q d<p 


V fi i 

^ 9x. t2 - ti 


+ 4 


8g* 


p 9 ip* 

+ ). i,, ~ — ; — r + i. = 0, 

V t-r~\ l+n 


8Z i( T 2 ) i,= i I ' 9z i < T 2> i+ P 


(24) 


and 


0 


y x * _i 

3=1 Z 3 9k i *2 “ ^ 


+ 4 


dip* 


-®SL_ + ?'i 

9z i ^ T i ) p =i P 9z i< T i> 


= 0 . 


(25) 


The last equation is unique in that it occurs at the left hand end point. 
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Just as we transformed the A* 
again to equation ( 25) giving p 


z i’ 


before we now do so 


9F 

9x. 

1 


+ An 


9g 

9x.(tj) 


V=1 


90 . 


v 9x^(ti) 


= 0. 


(26) 


Equations (24) and (26) are the first two equations from the trans- 
versality condition applied to zj. Equation (26) is in final form but (24) 
may be modified. This is done by combining it with the "next" equation - 
i. e. , since equation (24) is valid at the right hand end of the first arc, it 
would seem that we would be able to combine it with an equation which 
holds at the left end of the second arc (for the same x variable) giving an 
equation valid across the point t^. 


The next equation we need can be derived from equation (25) . The 
first point to be noted is that two of the terms in equation (24) vanish iden- 
tically. Earlier it was commented that, for our x variables, the range of 


definition was 


t. < t < 

J 


Vi- 


For our present transversality application we 


are dealing with the variables z 
these z's are defined for Tj ^ 


, . .v . for j = 1, . . . , m; i = 1, . . . ,n. Thus, 
(j-i)n+i * ’ 

t < r 2 and our variables are not defined at r 2 . 


So our g!' is a function of r 2 only at the point t (i. e. , right hand end point 
of the entire interval) . We may now write equation ( 24) as 


Z 


(3=1 


3_ 


k ~ h 


90 * 

9x. I 


2 

+ i. 

1+p 


= 0 . 


To obtain the equation which corresponds to the variable x., but on 
the next arc, we use a similar procedure. For the first term we have 


az* 

i+n 




00 * , ^ 80 * 
_£ dr ^ V ^ 


dt 


= 2 
13= 1 P 


9x. ta - tg 
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Bfa, 


for j = 2. So 


? * fa 

'? 3k l * - ** 


* 


+ *0 


3g* , ,* 3g* _ 

8 Vi (T2) ° 9z —- (Tl) 


n+i 


- £ 


dip* 


L * 


dip* 


Lj . v dz .(r 2 ) ^ v dz (tx) 

y=l n+i L v=l n+i' 1 


+ dz 


n+i 


i.. [z.(t 2 ) + Ax. . - z .(Ti)] = 0. 
l+p l i xj n+i 1 


Then for the point t* we have 


/* 


q 3$* a 

y x * _ £ _ i_ 

£ 3^ tj - t* 


+ AS 


3g* 


° 8 W T ‘' 


P 

+ T i 

A V 

v-i 


dip* 

7 - i. =0. 

8Z n+i (Tl) P 


Addition of this equation to the modified form of equation ( 24) cancels the 
constant <L . , , leaving: 

l+p 


*2 - *1 b ' X i 


z 


- z 


q A|_ 


_£ 


t3 - ta 3k. 


+ xs 


3g* 


8 Vi (T ‘ , 


p 

* 


dip* 


dz .(Ti) 
i>= 1 n+i 1 


= 0. 


Converting this to the original coordinates we have 


K fa 


fe” q 0 <* l t 2 + 

jki * 8 k i. 


+ An 


3g 


9x, 


Z* 


30, 


iW v 8x i (tz) 


= 0, 
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or 




8g 

3x.(t2) 


+ 


P dip 
y t — v__ 

tii v 8x i (t2, 



0. 


Now the point ^ is not unique; we can apply the same argument to 
any other point t. (j = 2, . . . , m-i) . Thus, 


9g 


W, 


0 9x.(t.) 


P 

9x. (t.) 
v=i i j 




9F 


9x. 

l 


= 0 


(28a) 


As a final datum from the transversality applied to the set of variables 
under discussion, we obtain a condition at the right hand end of the final arc. 
As was done above, we obtain 


9g 

9x.(t ) 
x m 


Z' 

v=i 


dip 

]v_ 

9x. (t ) 
i m 


9F 

9x. 


m 


= 0. 


(28b) 


Summary of Information Derived from the 
Transversality Condition 

Equations ( 22) , (23a), (23b), (26), (28a), and (28b) are the set of 
transversality conditions we have now derived for the problem. These are as 
follows: 


& + y i _ v ~ 

at i 


p 


n 


F - X x. ££- 
i=l 1 8x i 


= 0 


( 22 ) 


2 9^, 


x _2fiL + Vi _Ji + 

X ° 9t 8t 

m v=l m 


x . 9F 

F 'Z x i5T 

1=1 1 


m 


= 0 


(23a) 
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0 


(23b) 


-3S-+ Vi 
at, Lt ± v 9t . 


3 


v= i 



2 

i=l 


8F 


i 9x. 


t. + 

J 


t: 

3 


*0 


9g 

9x.(tt) 


+ 



v 9x.(tj) 



o 


(26) 


ag 


ax.(t.) 
i j 


p 

+ f. i 

i/=i y 


dip 

Y _v 

ax. (t.) 

i j 


9F 

ax. 



0 


(28a) 


8g 

8x.(t ) 
l m 


P 

2 ^ 

v=i 


dip 

9x.(t ) 
i m 


8F 

9k 



(28b) 


The conditions can be counted to determine the exact number derived. 
Only one each is obtained from equations (22) and (23a)', and, since j = 2,. . . , 
m-1, m-2 equations are derived from (23b). Thus, our first set yields m 
conditions. 

Equation (26) with i = 1, . . . ,n yields n conditions as does equation 
( 28b) . In equation ( 28a) , we have i = 1, . . . , n and j = 2, . . . , m-i, giving 
(m-2)n equalities. Our total from the second set is now nan, so we obtain 
a final total of m + mn = m(n+l) conditions. 

The next obvious question is the utilization of these equations. They 

are used for obtaining the values of the end point multipliers A. 0 , i , 

a total of pf i unknowns. Notice that this entire set of equations may be ^ 
written in matrix notation as 
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(where a summation is implied when a repeated index occurs) . Now since 
we have p+1 unknowns (with p < m(rH-i)) and m(n+l) equations we have 
more equations that unknowns (unless p = m(n+l)) and some equations had 
better be dependent; otherwise we have an inconsistency. That is, each 
(p4-2)x(p+2) submatrix must be singular; or, better, the augmented co- 
efficient matrix should have rank < p+2. 

(At this point Hunt argues that A 0 , A.^ do not simultaneously vanish 
and that Afi do not all vanish. His argument is primarily referenced to 
Bliss and will not be reproduced here since inclusion of the details would 
be a rather lengthy process.) 
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To proceed to a few definitions, we state first that an extremal for 
the problem of Bolza is an arc without corners in Rrj together with a set of 
multipliers zj(t) , ,A'£, (r) (with j = i,...,N;a = l,.».,M;0^r^ 1) 

such that the functions zj (r) , ( t) have continuous first order derivatives 

and satisfy the constraints qb >* = 0. An extremal is called non-singular if 
the determinant 


R* 


8 2 p» 

9z!8z* 
3 k 

8 <f> 

a 

8z’ 

k 


8 d>-\- 


3z! 

J 


0 


j » k - I, • • • )N 

a,a l9 = 1, . . . ,Q 


(30) 


is non-zero along the extremal. An extremal is normal if it satisfies the 
multiplier rule with a unique set of multipliers having A 0 = t. 


These definitions for the transformed problem can also be restated 
for the original problem we attacked. Thus, an extremal is defined to be 
an arc in R t and a set of multipliers, x^t), \ 0 , X^(t) (where x i (t) is 
the arc) with i = l,.,.,n;j = l m-1; /3 = 1 q; 


*3 


<t j+i 


have corners only at tj ( j = i, . . . ,m) . 


which 

(j) 


Between corners, Xj(t) and (t) 
have continuous first order derivatives, and the Xj’s satisfy the constraints 
^>(j) = o. An extremal is non- singular if the determinant 


8 2 f dcp y 



9x^85^ 9x^ 


| ij h — 1) • • • j u 

B -U 0 

8 k. 

1 


^*7 = ij • • ♦ » Q j 


(31) 


is non-zero along the extremal. An extremal is normal if it satisfies the 
multiplier rule with a unique set of multipliers having A 0 = 1. Normality 
for the original problem implies normality for the transformed problem, 
but this will not be proven here. Only the normal case will be considered 
in what is to follow. 



THE NECESSARY CONDITIONS OF WEIERSTRASS AND CLEBSCH 


The Weierstrass function, E, is defined by 


E* (t,z, z’,A* ,Z') = F* (t,z,Z',A* ) - F* (t,z,z',A* ) 


- 2 < z k - V Iff ■ 

k=i k 


(32) 


An arc is said to satisfy the Weierstrass necessary condition (for a minimum) 
if 


E* (t,z,z',\* ,Z') ^ 0. (33) 

at every element (r, z,z',A* ) of the arc for all sets (t, z, Z ') #(r, z,z') 
interior to the region Rf and satisfying the equations = 0. Every normal 
minimizing arc for the problem of Mayer or Bolza must necessarily satisfy 
this condition. 

Using ( 14) (which defines F* ) this condition may be written as 


E* (r,z, z', A* ,Z ') 


q 

-l)q+/3 [< ^(j-l)q+/3 


(t, x,X) 


(j-l)q+/3 


( t, x, x) ] 


(34) 


n 


^ lk 5f (Zk Z k )X (j-l)q+/3 


^ 9 < ^'(j-l)q+/3 ' 


9 \ 


t. . - t. 

J+i i 


46 


As has been done before, we indicate the derivative of a starred quantity with 
respect to the x^. This is done with the understanding that when we write 
< ^j-l)q+/3 (X ’* ,t) we <# > ^_ 1)q+/3 lx( T ,z),x(T,z,z') ,t( T >] , so that 

even though we can indicate the derivative of <p* with respect to x^, we 

really have the independent variable as t, not t. Thus, in <fts . and 

o. / . \ J ^ / tT p 




we substitute t = t. + (t. . - t.)r, x.(t) = z. . .. (r), 

j j-1 j 1 (j-l)n+i 


(l-l)n-H 

t . - t. 
J + i ~ J 


and t. = z. .. ... 

} (m-i)n+j 


To comment on the application of this test, we first note that we can 

consider all possible sets like (z,Z',t) in equation (32) , and then we pick 

out those which satisfy equation (32). Note that those we pick must satisfy 

the equations z' . , .= 0 and <£*.. .. „ = 0. Those we select are a sub- 

(m-l)n+j (j-i)q+i3 

set (of all possible choices). Those we vary to check the Weierstrass 

condition need only be members of the set of variables which correspond to 

the x's, namely z* . . .. We vary these one at a time; that is*, we choose 

(j-l)n+i 

one index where Z ' . .. ,, ^z.. .. .and choose Z* . .. ,.=z.. 

(j-l)n+i (j-i)n+i (j-i)n+i (j-l)n+i 

for all other values of j. Then the sets which satisfy equation (32) in z 

notation must satisfy equation (34) when we change to an x representation. 

In the other direction, if each of the m-1 expressions (j = 1, . . . ,m-l, 
one condition on each arc) is non-negative, then Weierstrass' condition is 
satisfied. 

We can transform these (m-1) equations to a complete representa- 
tion in our original notation to obtain the condition that 


Z^(t j+1 - t.) t,x,X) - ^ ]) (t,x,x)] 


“ L (X - x ) 

i=l 


8F(t,x,x, A) 


This equation must hold at every point (t, x, x,A) of E on the m-1 intervals 
t j < t < t. + ^ for all sets (t, x, X) ^ (t,x, x) which satisfy <p„ = 0 and are 
interior to our region R A . 
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By definition of our 
we must have 


function and the requirement that t ^ > t. 


E(t,x,x,A,X) = F(t,x, X,A.) - F(t,x,x,A.) 


n 


- 2 <*, - V 

i=i 


9F(t,x,x,A.) 

dx. 


S 0. 


( 35 ) 


This is the Weierstrass condition for the original problem and must hold 
for every minimizing arc. 

To delineate the Clebsch condition briefly, suppose we find 
numbers ir lt . . . , 7r n which are not all zero and which satisfy 

n d<p^\ t,x,x) 

A — Z-r. 7T. = 0 

H. dx. l 

1=1 l 


for all j = ,m-i. Then every normal minimizing arc must satisfy 

the condition that 


n n 


2 2 

i=l h=l 


9 2 F(t,x, x) 

9x.9x, 
l h 


7T 

i . 


> 0 


(36) 


as is shown in Bliss (Ref. 4) using the Weierstrass necessary condition. 


THE PROBLEM OF UNSPECIFIED DISCONTINUITIES 


We have thus far determined a set of four necessary conditions which 
treat the problem of an optimal trajectory with specified discontinuities at 
unspecified points. (Note, however, that we have not derived sufficient 
conditions. ) It is now possible to generalize the above formulation with little 
additional work to the point at which problems having unknown discontinuities 
at unknown time points can be treated. 



I 


In the case of a rocket vehicle which has constant thrust for each 
stage, prespecifying our comer points (i. e. , staging times) automatically 
prespecifies the mass that must be dropped. In order to determine the 
optimal weight staging we need a formulation which neither requires the 
staging times nor the drop weights to be specified. 

The difference may be illustrated as follows. For a fixed discon- 
tinuity we have an equation of the form 


x.(t+) - x.(t.“) = C, 
i J i J 

where C is a given constant. For a variable discontinuity we could write 

x.(t+) - x. (t. - ) = C(t. - t. ). 
i ) i J J J-l 

In the second case the discontinuity depends upon the length of the segment 
preceding the point t.. 

It is worthwhile to recall at this time that the discontinuity values 
assumed in the earlier part of this report, Ax^, did not enter into our 
boundary conditions. For this reason, equations (11) were specifically 
built into the transformed problem to account for the discontinuities. 

To proceed to the present problem, we now try to minimize 


J g[ t^, . . . , t^, x i( f i ) » • • • » x^(tj ), Xj( t^ ),..., 


X ( to ) , . . • , Xi (t ),..., X (t ) 1 

n * 3 ’ 1 m * n m 


(37) 

— §[ fi» • • • » ) » x (^2 )»•••* 1 

= gltj,. .. ,t m , x(t+), X(t£)] 
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where j = 1, . . . ,m-i and k = 2, . . . ,m. The cpp constraint equations re- 
main the same as they were previously. We can have three types of boundary 
conditions. The first, for no discontinuity is 


ip = x.(t.“) - x.(t. + ) = 0, (38) 

v 1 j 1 j 

which is stated at every corner where the x^'s are continuous. The second 
type of boundary condition is 


ip - x. (t. - ) - x. (t. + ) - Constant = 0, 
v 1 J 1 j 


(39) 


which is the type considered in the first part of this paper. Third, we have 


- x.(t. - ) - x. (t+) - D[ t u . . . ,t m , x(t+) , x( t~) ] = 0, 


(40) 


where the D's are specified functions (which, for example, give a relation- 
ship between the burning time of a stage and the drop weights) . 


The basic difference between the present formulation and the previous 
work was that we had to include the constraints (11) which related to the dis- 
continuity Ax^. If we explicitly include the constraints given above by equa- 
tions (38) , (39) , and (40) , then we can exclude the constraints of equations 
( 11) . (What we have really done is to convert all constraints to the form of 
equations (11) and then renumber them from the beginning as ip p , 
v = 1 p^m(n+l) . 


The modification due to excluding the constraints of equations (11) is 
rather minor. In the development of equation (24) , an end point multiplier 
i. turned up which was then eliminated by the same term in a subsequent 

equation. This time we do not get the term f to occur in either case and 

there is no necessity to add the equations in question to eliminate it. 

Notice that we have both the points tj + and tj“ in the constraints ip v since 
our discontinuity is not fixed. Previously it was argued that the terms of 




the form \q 


9g 


9z.(t 2 ) 


and 




occurred only at the point t m . In the 


present case this is not true since the discontinuity is not fixed. 



vVlthout further ado, we write the two equations that replace equa- 
tion (28a) as 


Q P dib 

\ 8g , , V e JL 

0 ax.(t + 7 £j z v 8x.(t+) 
i J v=l i j 


8F 

ax. 


V 


= 0 


*0 


M. 


p 

+ fjJL 

u . v 




dx.(t.-) Lj , v dx. (tr) 
1 J v=i 1 J 


8F 

ax. 


]t f 


= 0 


The D functions occurring in equation (40) modify the problem by their 
occurrence in the ip constraints. 


RESULTS AND CONCLUSIONS 


Our first set of variables were x< 


,x . These are the coordinates 


— > - n - 

of the point in phase space. After transformation these become 

^1> • • • J J \ 4 

1 (m-i)n 

The second set of variables were fy, . . . , t^. These are the "times" 
(i. e. , values of the independent variables) at which discontinuities occur. 


After transformation these become z 


these two, we now have variables 
tion, and we set ( m-1) n+m = N. 


(m-i) n+1* 
i» • • • » z 




Between 


(m-i) n+m 


(m-i) n+m* 
after the transforma- 


The first set of constraints we labeled as 


(j) _ 


= 0, where the j 


indicated the particular "arc" (or interval between two values of tj) on 

which the constraint applies and the /3 indicates the particular constraint 

on this arc. After transformation these become d>t* . . . , <£*> , , (q<n) . 

1 (m-l)q 


The second set of constraints were z* .. . = 0, which were intro- 

( m- i) n+j 

duced to account for the time points being constant. Our total set 

of constraints were <p% ... .dr] .. , and we set (m-i) q+m = Q. 

1 ( m- 1) q+m 


The first set of end conditions were the ip v = 0 constraints, or 
ipi ( p=£m(n+i) ) . After transformation these became #*»...,$ p . 
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The second set of end conditions accounts for the discontinuities of 

the Xj's at the ti points which are, after transformation, z,. . > ,.(1) -Ax. .- 
i J (j-l> rn-i ij 

z.. ,, . = 0, which were written as ib*. . . =0. 

(j-l) n+n+x r (j-l)n+i+p 

The third set were t, ? # ( m _ 2) wpf j = ° and r 2 - 1 - * (m . 2) ^ = 0. 
This gives the total set as $£, . . . » ^( m _2) n+p+i’ an d we lab®! (m-2) nfp+l = P. 

Next, certain Lagrange multipliers S. pt A. 0 were introduced and 
the function Ps< formed. Application of the multiplier rule yielded the 
following results: 

The Euler equations for the original set of x variables are 


¥~ ~ f , = 0 (t. t < t ). 

3x, J 3x, (j-i)n+i J j-1 


1 t. i 


Transversality conditions: 


9 1 


+ E \itf- 

1 v=i 1 


n 


9F 


V • oi 

- / , X , TT 

H a 1 9x. 
1=1 


= 0 






m v=i 


n 


F - 2 x. ££■ 
h. i ax. 
1=1 i 


m 


= 0 


Q P dip 

, 9g_ . y i Jjl _ 

x ° at. LtK 9t . 

3 V=i j 


y • sf 

H, X i 8x. 
1=1 


t? 


= 0 


3F 


Q p dip 

— 8 .g. + >' a —JL. 

9x i(ti) v 9x.(t 1 ) ax. 


= o 
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*0 


*0 


9g 

y ^ 

Q? 

9F 

* 

dx^ tj 

1 v 
v=t 

9x.(t.) 

CO 

*f 

9g 

P 

y § 

dip 

Y v 

8F 

*m 

9x.(t ) 
l m 

U A V 

8x. (t ) 
l m 

8k. 

l 



= 0 


= 0 


Weierstrass condition: 


n 


F(t, x,X,A) - F(t,x,x,\) - Jj (X, - x ) 

i=i 


. , 8F(t,x,x,\) 


8x. 


s 0. 


Clebsch condition: 


V V 9 2 F (t, x, x) 

H. t*. dk.dk, 
i=i h=i l n 


7T.7T, — D. 
i h 


The problem with unspecified discontinuities is similar to equation 
(40) in which the D functions appear, modifying the ip p constraints. Aside 
from this difference, the only modifications occur in the fourth equation 
listed above under transversality. This becomes the two equations 


p dip 

-§s._- + y i — il - 

8x7(t*j v 9x. (t. + ) 

1 ] v=i 1 ] 


8F 
1 8k. 


"j 


= 0, and 




dg 


P 

+ 2j K 


8 1 


8 *i<y> ’ & v dx iV 


9F 

8x. 


"i 


= 0 


The Weierstrass and Clebsch conditions for this case have not been 
investigated here. 
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